Super-slippery Carbon Nanotubes: Symmetry Breaking breaks friction 
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The friction between the walls of multi-wall carbon nanotubes is shown to be extremely low in 
general, with important details related to the specific choice of the walls. This is governed by a 
simple expression revealing that the phenomenon is a profound consequence of the specific symmetry 
breaking: super-slippery sliding of the incommensural ^ " ' ' ^^■^■'1 11 

principles of tribology, offering a recipe for the lubric; 



PACS numbers: 61.48. -l-c, 62.20.Qp, 64.70.Rh 



Various possible applications of the carbon nan- 
otubes [0 motivated extensive investigations of their re- 
markable physical properties, from the mechanical to the 
optoelecronic ones. The symmetry based prediction |^ 
of the extremely small friction between the walls of a 
double-wall tube (DWT) has been confirmed recently 
both by the thorough numerical analysis ||^ and the 
experiment [Q. Here we give complete theoretical dis- 
cussion of such a behavior. To this end we analyze 
the DWTs W'-W (Tab. |). The waU W is (12,12); 
its radius is D'/2 = 8.14A. The waU W is one of the 
remaining tubes (711,712). The first three with radii 
Z5/2 « i?'" ~ 4.7A, suite as the inner wall, and the last 
four with D/2 « = 11.57A, as the outer one. The 
first part of the paper gives the exhaustive insight to 
the variety of the interesting phenomena. Then the re- 
sults are reconsidered to extract three general principles 
profoundly relating interaction and symmetry, leading to 
the final conclusion: the reduction of friction is due to 
the incompatibility of the DWT walls symmetries. 

Any single- wall tube (ni, 77.2) is completely determined 
by its symmetry group The [/-axis maps the atom 
Cooo to the second one Cqoi- Then the successive ro- 
tations C„ for the tube (p-period (j) = 27r/77 (71 = 
5(771,772); Q - the greatest common divisor) yield the 
initial monomer with 27i atoms Cqsm (s = 0, ... ,77 — 1; 
7t = 0, 1). Finally, the screw axis generator (C^jTja/g) ar- 
ranges helically the monomers, translating (for na/q) and 
rotating (through 2Trr/q) them along the z-axis. Here, 
q is the order of the isogonal group principle axis, a is 
the translational z-period and r is the integer related to 
the chirality. Let the x-axis of the coordinate system is 
to be rotated for $ and translated for Z along z-axis 
to coincide with the tube J7-axis. Then the coordinates 
I'tsu — {D/2, (ptsu, Ztsu) of the i-th monomer atom are: 

Vtsu = (-l)"v3o + 27r(- + -) + $, (la) 
q n 

Ztsu = {-iTzo+t-a + Z, t = 0,±l,.... (lb) 
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Here ao = 2.46lA, 

Firstly, we look for the potential V{lp, z) produced by 
W' at the arbitrary point r = (_D/2,<p, z). Since neces- 



FIG. 1: The potential V{<p, z) of (12,12) at p = 7?'" (left) and 
p = (right) within the ip- and z-periods <j)' and a' , with 
the expansion amplitudes below. 



sarily invariant under the symmetry group of W, V{ip, z) 
can be expanded in the basis of invariant functions, har- 
monics. The suitable basis invariant under the f/-axis 
defined by $ and Z is [M = 0, 1, . . . ; w real): 

C^\lP, z) = cos{M{ip - $) + 27rw(z - Z)), (2a) 

There is unique combination of C^' and C^^^, namely 

A'J{^,z) = C^'{^,z) + C^l{^,z), {M,uo>Q), (2b) 

invariant under the mirror planes through U . The roto- 
helical symmetries impose restrictions on M and uj; in- 
deed, the invariance under C„ and (C^jna/g) implies: 

7\/ = mod 77, (rotations), (3a) 
Mr + ujna — Q mod g (helical). (3b) 

The solutions of this and the other systems appearing 
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TABLE I: Family of DWTs (m , n2)-(12, 12). The row of W= 
(ni,n2) gives: symmetry group of W (line group Tg{a)Dn 
or T^{a)D^h for W chiral or achiral) and of W-W (Gn), 
common solutions X^' of ^ and ^ in the form of (^. For 
commensurate DWT Gn = Tq(A)Cn (line group), Si — j = 
q/a, S2 = Y = Q / A and S = 77'/r. For incommensurate 
DWT Gn = Cn (point group), = q, S2 = G{q,q') and 
S — qq'Q{q, q ). The translational periods a and A are in the 
units of ao- 



Tube Line Group Gn 


X^i' S*! S 


(12,0) T^i{V^)Dr2h C12 
(7,7) TUl)DrH T/(l)Gi 
(11,2) Tii{7)Di T/(7)Gi 


ig* 24 24 24 
if** 14 1 336 
7f** 14 i 2352 


(12,12) T2^4(l)-Di2h 


24 


(17,17) Ti^{l)Di7h T/(l)Gi 
(18,16) Tle^ V65T)D2 G2 
(29,1) Ti%' (V2613)£>i Ci 
(24,9) TiM{V29l)D3 G3 


if** 34 1 816 
868 4 5208 
^20904 -^7^2 2 20904 
if^* 582 6 2328 



latter on will be given in the form (L, K = 0, ±1, . . . ) 

M^M*L, Lu^uj*K, YL + K = mod X, (4) 

with specified M* , uj* , X and Y (e(X,r) = 1). In the 
rectangular 2D lattice with the periods M* and oj* , the 
allowed pairs (Af , a;) form the sublattice with the cell 
(0,0), {0,-uj*X), {M*,-uj*Y) and {M*,uj*X - uj*Y). 
This cell with area M*uj*X contains only trivial solution 
M = Lu = 0. Particularly, (|^) is solved by (^) for uj* ~ 
a^^, M* = n, X = q/n and Y = r; the cell area is 7 
q/a. This singles out the harmonics C^^ and A^^ of the 
chiral and achiral tubes. The harmonics (p- and z-periods 
arc 4>/L and a/lXl; those with M ^ n and to = ±a~^ are 
invariant under the roto-helical symmetries of W only, 
while the others have L\K\ times finer periodicity. 
Consequently, the allowed M and oj for (12,12) tube. 



M = 12L, 



K/ao, L + K = mod 2, (5) 



single out the harmonics ^'^ in (pbl). Further, for the 
W wall we fix $ = Z = (x-axis is chosen as the U -axis 
of W). Therefore, one finds: 

Vi^,z) = J2Vc{r~rU- J2' a'jA'^'i^^z), (6) 

tsu M,u)>Q 

where the prime indicates the summation over the so- 
lutions of (P) only. We calculate the potential by the 
left part of^6[) with Vc being the Van der Waals type 
interatomic potential fitted to the 7r-'--bonding The 
sum over 82 monomers (41 elementary cells with 1968 
atoms) is scanned in 41 x 41 points within ip G [0,7r/6) 
and z € [0, 1). The fast Fourier transform yields the am- 
plitudes a^tx^^ {L,K ^0,..., 20) on the right of (|). 
The results for i?''' and are in Fig. |l|. As the po- 
tential V{ip, z) practically does not vary on 10~^A scale 



(in the vicinity of i?™ and these plots refer to all 

the DWTs listed in Tab. |. Arbitrary constant ag is set 
to zero. No aliasing occurs, the functions excluded by 
(^) are with zero amplitudes, while the amplitudes of the 
harmonics rapidly decrease with L and K, and become 
negligible (within the numerical error of lO^^^meV) out 
of the depicted range. This property is nicely illustrated 
in the density plot (Fig. ^ of the potential V{lp, z): It is 
notably periodic only with the periods of the tube. 

With fixed position of W'-wall, the W'-W relative po- 
sition and mutual interaction are determined by $ and 
Z of the W-wall [/-axis. The interaction is the sum of 
the potential (^) values in the positions (|^) of W atoms: 



A/,w>0 tsu 



M a'M 



{'■Ptsu, Ztsu)- 



(7) 



For the infinite W-wall the summation over u, s and t 
easily gives the interaction per the W-wall atom: 



Vooi'^,Z) = E" cos(A/$ + 27rcjZ) X 

M>0,ui 

X 2a|^j| cos(A/i^o + 27rwzo)- (8) 



Here the double prime restricts the summation to M and 
Lo satisfying simultaneously (^) and ^ for the W-wall 
parameters. This means that the interaction is mediated 
only by those W'-harmonics which are also invariant un- 
der the roto-helical symmetries of W. Then the allowed 
M-LU values (solutions of the system (||)-(|5|)) are given 
by (|) with y = 1 and M*, lu* and X hsted in Tab. |. 
The absolute value of the constants in the second line of 
(^ roughly determines variation of v{^,Z). Thus, the 
corrugation |^ over W-atom is small whenever the area 
A!*uj*X is large: large M and lu indicate small \a^^\. 

For the incommensurate (a'/a irrational) walls lu* = 
0. Consequently the mediating harmonics are constant 
along Z, yielding Z-indepcndent Voc- Further, M* = 
qq' /G{q, q') {X ~ 1 imposes no restriction). The ideal 
infinite walls slide along the tube axis exactly without 
friction, while the rotational friction oscillates with the 
<i>-period 2-k/M* . Fig. ^ presents the calculated interac- 
tion Woo(*i>) for the (12,0) wall; its ^-period is (j)' /2. In all 
other cases M* is rather large, allowing only very high M 
of the mediators; their amplitudes are beyond the acces- 
sible numerical precision and the computed interaction 
vanishes, indicating small friction. In the commensurate 
cases M* is large: the ^-dependence of the interaction 
stems from the harmonics with extremely small ampli- 
tudes, giving negligible rotational friction. The calcu- 
lated interaction is only Z dependent (Fig. |^), with the 
Z-period a'/2 (since w* = 2). Particularly, for the (11,2) 
wall X = 7. This additionally forbids all the numeri- 
cally significant harmonics: (11,2)-(12,12) walls friction 
is extremely weak. 

For the finite W-wall with m monomers, the interac- 
tion is also obtained from (^. The horizontal symmetry 
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axis of W is in the middle of the ring: $ 

r-7 rz , n m — 1 mi i ■ 



is multiple of the W-W translational period A, all the 
numerators vanish. Oijy-Ihe infinite wall mediators con- 
S akc infinite wall corrugation/atom 
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FIG. 2: The corrugation/atom dependence of the VK-wall 
monomer number for the considered DWTs. Inset: Woo($) 
for the wall (12,0), and «oo(^) for (7,7) and (17,17) (for the 
other walls t;oo(4>, Z) is below the precision). 



Analysis of (y) reveals quite interesting phenomena. 
The second line in (jp) apparently suggests 1/m dump- 
ing of the corrugation/atom quasi-oscillations. However, 
both the numerator and denominator vanish for the in- 
finite walls mediators (H), giving TO-independent terms 
which are summed in Vao- The additional exclusively 
finite wall mediators give the dumped oscillations su- 
perposed to foo, as transparent from Fig. |^. Although 
numerically v^o may vanish, really it is only small, and 
the total corrugation always increases with m. For the 
(7,7) wall only the infinite wall mediators are numerically 
recorded, while for (17,17) they prevail by several orders 
of magnitude. Since they are with M = 0, the rotational 
friction is almost negligible for any m. On the contrary, 
for (12,0) the rotational friction is dominant. For the 
other tubes the infinite wall mediators are not recorded, 
making both rotational and translational friction compo- 
nents significantly oscillating with m and approaching to 
zero. Peculiarity of the commensurate cases is that when- 
ever TO is multiple of qA/na, i.e. when the ring length 
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T>f V \mnrianT imaer tm sfmmm-y vj Vrr Tilus, only 
M n ITtSl^S'^ction of V to the subspace of the W^- harmonics 
counts, while the contribution of the orthogonal part van- 
ishes. Indeed, the interaction is mediated only by those 
W' harmonics that are additionally invariant under the 
roto-helical symmetries of the other wall. Even for them, 
from (|) follows that of W is orthogonal to all C^',' 
of W unless M — M' and lo ~ . Therefore the pro- 
jected is Cf„((p, z) cos(M$ ± uZ). Due to the W- 
wall invariance of C*^, its value is the same C^'^((^Oj2o) 
over all the W atoms Gtsu- This completely explains (||). 

Taken together, A and C show that the symmetries of 
the interacting systems W and W restrict their inter- 
action only to the W'-harmonics having invariant parts 
with respect to the W-wall symmetry. This is performed 
in two consecutive steps. Firstly, the roto-helical symme- 
tries of the W and W walls are collected into the sub- 
groups (i'^^)-type [|j) G and G' of mutually commuting 
elements. This enables the choice of the harmonics of 
G' being either harmonics of G (becoming mediators) or 
orthogonal to them (giving no contribution). Secondly, 
as the parities (J7-axes) do not commute, they select (for 
each of the symmetry groups) the linear combinations 
(^ of G'- and G- harmonics, which are not mutually or- 
thogonal. The parities thus influence the interaction by 
the cosine terms which dump the mediators amplitudes 
in (^). Analogously, the first line in (|^) is also the dump- 
ing cosine, while the second one is the sum of C^^ values 
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at different monomers (wlien the ring symmetry group 
Dn acts on a single atom, a single monomer is obtained; 
consequently, the £)„-harmonic has constant value 
only along a single monomer). 

Hence, only the mutually commuting symmetries pro- 
duce the rarefaction of the mediators (explicated by the 
double primes in (H) and (^)). The corrugation is mainly 
determined by the phenomenological harmonic ampli- 
tudes (reflecting the nature and intensity of the inter- 
atomic forces), i.e. by the commutative symmetries of 
the walls, selecting the corresponding mediators. Still, all 
the symmetries contribute to the fine tuning of the inter- 
action: the dumping factors and the mediators values at 
the atoms result in the interaction position dependence. 
Just this causes the friction. 

Clearly, the rarefaction is the most important effect. 
Its estimation in terms of the commuting groups quite 
profoundly clarifies the interplay between symmetry and 
friction. According to the note after (^), the cell area of 
the wall harmonic sublattice is 7 = q/a. This quantity 
directly measures the roto-helical symmetry: the more 
symmetric the wall is, the severer are the restrictions 
(^, and the harmonics are sparser with the larger cell. 
(Alternatively, since all W-atoms are obtained by the ac- 
tion of its symmetry group on Cqooj their linear z-density 
measures the order of the group; and 7 is the half of this 
density.) Further, the area S = AI*ui*X points to the 
lowest mediators, giving the estimate of the inverse fric- 
tion. Simple inspection of Tab. || in the commensurate 
cases shows that S — 77Vr, where T = Q/A measures 
the DWT symmetry group. The established result in the 
form S = \G ® G'\/\Gr.\ (the order of the direct prod- 
uct is \G ®G'\ = |G||G'|) means that S is the symme- 
try breaking rate from the group G ® G' to the actual 
DWT symmetry group Gn = G n G'. And, G ® G' 
is the symmetry group of the non- interacting systems: 
independently performed transformations of G and G' 
preserve the system energy. 

As we have shown, the interaction itself imposes a spe- 
cific symmetry breaking. In the usual manner this pro- 
duces S = |G||G'|/|Gn| equivalent minima in the to- 
tal energy (recall Jahn- Teller-effect, crystal phase transi- 
tions, particle physics symmetry breaking). While the 
minima density increases with the breaking rate, the 
variation of energy between them, proportional to the 
mediators amplitudes, decreases according to the prin- 
ciple B. In the extreme cases the energy is constant 
along some paths. This Goldstone super-slippery mode 
of the incommensurate walls is the phason of the incom- 
mensurate phase transitions. Even then the remaining 
pure rotational friction reflects the breaking of the rota- 
tional symmetry. The z-independent mutual interaction 
iuj = 0) makes the walls effective symmetry groups Cq' 
and Cq (isogonal to G' and G), with the intersection 
^G{q,q'); ^^^^ breaking rate qq' /G{q,q')- For the incom- 
mensurate DWT this equals M* for (||) (Tab. |), giv- 



ing the first allowed mediator (M*,0). With M* = 24 
for (12,0)-(12,12) the walls therefore interact through the 
W'-harmonic with L = AI/n^= 2. Theoretically, this is 
the lowest possible DWT mediator, yielding the maximal 
rotational W-W' friction. To conclude: in the view of 
the one-to-one correspondence of the single-wall tubes to 
their symmetry groups, the symmetries of the walls are to 
the large extent incompatible; the tremendous symmetry 
breaking and harmonics rarefaction leaves the mediating 
harmonics with very small amplitudes. 

The performed analysis and conclusions are essentially 
general, and refer also to other structures with different 
type of symmetries For example, for DWT with both 
walls incomplete, the corrugation should decrease with 
the rotational symmetry breaking rate nn' /N. While A 
and C are purely group theoretical (the proof is based on 
the group averages), B is phenomenological. Its elab- 
oration is beyond the scope of this letter. Note that 
only its global validity is expected: the amplitude of the 



harmonics may locally increase (e.g. \a' 
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for 



D/2 ~ i?°"*), revealing the complex structure of the sys- 
tem or some hidden symmetry of the potential. 

Several remarks we add at the end. The minima of 
the potential (§) determine the equilibrium positions 
of the W-wall. For the achiral walls minima of Voo 
(Inset of Fig. |2|) show that such DWTs have mirror 
planes Linear and angular momentum quantum 

numbers of the vibrational modes of the rotational and 
translational walls displacements are /c = to = 0; such 
modes can be detected by Raman scattering. Their fre- 
quencies may be estimated expanding the potential at 
minima For the considered sufficiently long DWTs 
the maximal ones are given by the Inset of Fig. |^: 
the rotational lj^^^''^'' « 21cm^^, and two translational 

cj^''^'' ~ ISlcm^-'^, cj^^'"'^'^'' « 85cm^-'^. The others (al- 
most) frictionless degrees appear as the (almost) acous- 
tic modes. This may be a hint to test the friction by the 
Brillouin scattering. Combined with the described reso- 
nant corrugation effect, this can facilitate the experimen- 
tal identification of the walls of the sample tubes. Finally, 
the rotational and translational displacements are in gen- 
eral coupled in two helical normal modes, corresponding 
to the maximal and minimal gradient of the mediators. 
This interesting (and possibly applicable) screwing effect 
may be estimated within the presented framework. 
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